Specific Heat of the S — 1/2 Two-Dimensional Shastry-Sutherland Antiferromagnet 

SrCu 2 (B0 3 ) 2 in High Magnetic Fields 



o 
o 

(N 
S3 



c3 



I 

S3 
O 

o 



> 

On 
O 

in 

o 
m 
O 

03 



S3 
O 

o 



X 
S3 



H. Tsujii, 1 R. C. Rotundu, 1 B. Andraka, 1 Y. Takano, 1 H. Kageyama, 2 and Y. Ueda 2 

'Department of Physics, University of Florida, P. 0. Box 118440, Gainesville, Florida 3261 1-8440 
2 Institute for Solid State Physics, University of Tokyo, Kashiwa, Chiba 277-8581, Japan 

(Dated: February 6, 2008) 

We characterize the field-induced magnetic phases of SrCu2(B03)2, a frustrated spin-1/2 Heisen- 
berg antiferromagnet in the two-dimensional Shastry-Sutherland lattice, using specific heat in 
magnetic fields up to 33 T. We find that the spin gap persists above the expected critical field 
H c — A/g/iB of 21 T despite the appearance of magnetic moment in the ground state. At the 
magnetization plateau at 1/8 of the saturation, the S z = +1 triplets that carry the magnetization 
of the ground state are observed to form a two-dimensional spin gas of massive bosons. A spin gas 
consisting of the same number of massive particles continues to completely dominate the specific 
heat in the field region above the plateau, although the magnetization increases with increasing 
field. Ordering is observed at a temperature immediately below the spin-gas regime. 

PACS numbers: 75.45.+j, 75.40.Cx, 05.30.Jp, 67.40.Db, 75.50.Ee 



Geometrical frustration is one of the important sub- 
jects in strongly correlated systems pj. In particular, 
the fascinating interplay of geometrical frustration and 
quantum fluctuation in two dimensions provides a fertile 
ground for new physics, as has been theoretically demon- 
strated for Heisenberg antiferromagnets in the triangular 
lattice and the kagome lattice. For the 5 = 1/2 triangu- 
lar antiferromagnet, there now appears to be a consensus 
that Neel order must exist despite the strong frustration 
0. The S = 1/2 ka gome antiferromagnet has been pre- 
dicted to possess a gap separating the ground state from 
upper triplet levels and a band of low-lying singlet exci- 
tations within the triplet gap, although the exact nature 
of the ground state is still under debate Q . These pre- 
dictions are yet to be borne out by experiment, because 
no isotropic S = 1/2 model system has been identified in 
the laboratory for these lattices. 

In the last few years, there has been a rapid progress 
in the understanding of the S — 1/2 Heisenberg antifer- 
romagnet in the Shastry-Sutherland lattice Q, another 
two-dimensional frustrated geometry, for which the exact 
ground state is known at zero and low magnetic fields de- 
spite geometrical frustration. The progress owes largely 
to the discovery of the spin gap and other novel mag- 
netic properties in SrCu2 (1303)2 0, the only known lab- 
oratory model for the geometrically frustrated S = 1/2 
two-dimensional Heisenberg antiferromagnet. 

The low-field ground state of the Shastry-Sutherland 
antiferromagnet is simply a product of spin-dimer sin- 
glet wave functions, and the dimerization by the diag- 
onal bond J causes a gap to form 4]. The gapped 
triplet excitations, which are also spin dimers at least 
to a good approximation, have a large mass as a result of 
the unique lattice topology involving orthogonal arrange- 
ment of dimers I n SrCu2(B03)2, the coupling 
strengths extracted from experiments are J = 85 K for 
the intra-dimer exchange and J' = 54 K for the inter- 
dimer exchange 10] or J = 70 K and J 1 = 42 K pj|. The 
large J' makes the energy gap A substantially smaller 



than J 0, [^. In fact, the ratio J' /J is very close to 
the critical value 0.70 |E IE El for the quantum phase 
transition to a square-lattice Neel state. In other words, 
J and J' are extremely frustrated. 

The striking feature of SrCu2 (BOa)2 is the existence 
of magnetization plateaus at 1/8, 1/4, and 1/3 of the 
saturation magnetization 0, Although there is a 

substantial spin-phonon coupling in this compound |15| . 
the plateaus arc an intrinsic property of Heisenberg spins 
in the Shastry-Sutherland lattice and due to geometri- 
cal frustration [(J, Il6| . The spin state at any magnetiza- 
tion plateau is obviously gapped , and a possible con- 
nection between magnetization plateaus and quantized 
Hall-conductance plateaus has been pointed out 0. In 
contrast, the states between H c and the first plateau, 
between two adjacent plateaus, and between the high- 
est plateau and the saturation field have been believed 
to be all gapless by a prevalent view on magnetization 
plateaus. 

In this paper, we examine three field-induced magnetic 
phases of SrCu2(B03)2: the one between H c and the 
first plateau with 1/8 of the saturation magnetization, 
the plateau phase, and much of the phase between the 
1/8 and 1/4 plateaus. The results refute the conven- 
tional view on the nature of the first phase and uncover 
the presence of a two-dimensional spin gas in the other 
phases, as well as ordering in the third phase. 

Prior to our work, two-dimensional Bose gas has been 
realized in spin-polarized atomic hydrogen adsorbed on 
the surface of superfluid 4 He [l^. However, the spin gas 
in SrCu2(B03)2 is the first example that is amenable to 
thermodynamic studies. 

The calorimeter for this investigation employed the re- 
laxation method [l9| . In order to optimize the relax- 
ation time at different temperature and magnetic-field 
regions, three samples of different sizes were used, all 
taken from a high-purity single crystal whose growth pro- 
cedure has been described in Ref . [2(| . The sample sizes 
were roughly 3 x 3 x 0.07 mm 3 , 4 x 3 x 0.3 mm 3 , and 
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FIG. 1: Specific heat of SrCu 2 (B0 3 )2 at 14 T, which is below 
the expected critical field H c of 21 T, and at 22 T and 24 T 
in the region between H c and the first magnetization plateau 
at 1/8 of the saturation. The solid lines are the fits discussed 
in the text. The broken line is the T 3 phonon contribution 
that has been subtracted from the data. 



3 x 3 x 1.3 mm 3 , and the magnetic field was applied par- 
allel to the crystalline c axis, perpendicular to the two- 
dimensional layers of the S = 1/2 Cu 2+ spins. In this 
field orientation, the spin gap for the S z = +1 tri plet s is 
split into two branches that are 4.1 K apart J 2ll |22 | by 
the Dzyaloshinskii-Moriya interaction 123 . [24J . Accord- 
ing to the g factor measured by ESR [21j , the lower gap 
is expected to close at H c = A/g/i^ of 21 T. 

The measurements were done primarily in magnetic 
fields between 22 T and 33 T at temperatures between 
0.53 K and 15 K. Additional data were taken at zero field 
and 14 T for comparison with earlier specific-heat results 
at magnetic fields up to 12 T [l^. At zero field, there 
was excellent agreement between the present data and 
the earlier result, indicating a high sample quality and 
reproducibility. The 14 T data precisely follow the field 
dependence observed in the earlier experiment, with the 
spin gap decreasing linearly with magnetic field and the 
position of the broad maximum remaining at 8 K. 

Figure 1 shows the specific heat at 22 T and 24 T, 
where the spin gap is expected to be absent and the 
ground state is magnetic [5J, [lj] , together with the 14 T 
result for comparison. The T 3 phonon contribution to 
the specific heat has been subtracted from all the data 
presented in this paper, as determined by Kageyama et 
al. |l2j| Instead of a broad maximum at 8 K, a new 
hump appears at 22 T and 24 T at temperatures below 
1.4 K. The position of this low-temperature maximum 
decreases with increasing magnetic field. At 22 T, the 
low-temperature tail of the hump clearly shows a gapped 
behavior, whereas the 24 T data do not extend to a suf- 
ficiently low temperature. 

The overall behavior of the specific heat can be ex- 



plained in terms of the dimer-quadrumer model of Hof- 
mann et al. |25| , who incorporated the repulsion between 
triplet spin dimers and bound pairs of triplet dimers as 
a mean field. The line for 14 T in the figure is a calcula- 
tion using exactly the same parameters chosen by these 
authors to fit the specific-heat data at fields up to 12 T. 
Similar to their fits, the 14 T calculation reproduces the 
experimental data very well at temperatures below the 
broad peak. 

In order to fit the data at 22 T and 24 T, we have 
found that it is necessary to allow the gap energy A + of 
the S z — +1 triplets and the mean- field strength vo of the 
repulsion term to be fitting parameters [2(| . The best fits 
shown in the figure are obtained by choosing A + = 4.6 K 
and 3.2 K for 22 T and 24 T, respectively, and «o = 
10 K for both fields. The fit at 22 T is good primarily 
at temperatures below the broad low-temperature peak, 
whereas the fit at 24 T agrees well with the data up to 
4 K. The good overall agreement leads us to conclude 
that the ground state is gapped at both fields. 

The remarkable persistence of the gap in the field- 
induced magnetic phase contradicts the widely held view 
that the ground state above H c evolves with increas- 
ing magnetic field by simply increasing the number of 
S z = +1 triplets within itself. We propose in contrast 
that mixing of at least one of the S z — +1 triplet levels 
with the singlet ground state is responsible for the evolu- 
tion of the magnetization below the 1 /8 plateau and for 
the anticrossing that keeps the gap open even above H c . 
The mixing is probably caused by the Dzyaloshinskii- 
Moriya interaction. The gap keeps the ground state in 
a spin-liquid phase with a non-zero magnetization. This 
explains why the translational symmetry is maintained 
in this field region even at a temperature as low as 35 
mK, as has been found by NMR [27|. We conclude that 
the region from zero field up to the lower edge of the 
1/8 magnetization plateau is a contiguous phase, with 
no ordering even at fields above H c . 

The situation resembles the behavior of the Haldane- 
gap antiferromagnet NENP, where the staggered g tensor 
of two non-equivalent Ni sites causes the triplet excitation 
to anticross the ground state at H c [2^, |2jj and magneti- 
zation develops with no long-range order due to the inter- 
chain coupling 30] . It is interesting to note that isolated 
triplet spin dimers in the Shastry-Sutherland lattice have 
been predicted to be unstable below the 1/8 plateau and 
they cannot be constituents of the ground state in this 
field region |3l| . 

The present result is consistent with the ESR experi- 
ments [21], [23] , which have found no excitation level that 
reaches the ground state at H c = 21 T. In fact no such 
level has yet been discovered up to 40 T. Instead, the 
behavior of the lower S z = +1 branch, and possibly also 
the upper branch, suggests anticrossing around H c [22| . 

The fitting parameters are also in good agreement with 
the ESR data H2, which give A+ = 4.8 K and 6.6 K for 
the two branches of excitations at 22 T. The ESR values 
at 24 T are 5.2 K and 5.5 K. Our value for the repulsion 
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FIG. 2: Specific heat of SrCu 2 (B0 3 )2 at 27.5 T (o), which is 
the midpoint of the 1/8 magnetization plateau, and at higher 
fields in the region between the 1/8 and 1/4 plateaus: 29.6 T 
(□), 31 T (V), 32 T (o), and 33 T (a). The dotted line 
indicates C = R/8. The solid line through the 27.5 T data 
points below 5 K is a single-parameter fit for a non-interacting 
two-dimensional Bose gas. Other solid lines are guides to 
the eye. The T 3 phonon contribution shown by the broken 
line has been subtracted from the data. The inset shows the 
ordering temperature as a function of magnetic field. 



parameter vq — 10 K at 22 T and 24 T is considerably 
smaller than 17 K for fields up to 12 T 25] and at 14 
T. This is expected, since the mixing of an S z = +1 
excited state with the ground state at fields above H c 
turns it into a more singlet-like object, which experiences 
less repulsion due to the inter-dimer exchange J'. 

Figure 2 shows the specific heat for fields of 27.5 T and 
higher. 27.5 T corresponds to the midpoint of the 1/8 
magnetization plateau, where Cu NMR has found evi- 
dence for a rhomboid ordered structure involving S z = 
+1 triplets 27j. Within our experimental accuracy, the 
specific heat at this field is 1/8 of the gas constant R 
over a wide temperature range except for a gradual rise 
at temperatures above 5 K and a small decrease at tem- 
peratures below 1 K. This is conclusive evidence that 
a two-dimensional spin gas is formed by the S z = +1 
triplets whose number equals 1/8 of the total spin dimers. 
This conclusion is consistent with the n B NMR spectra 
|27|. according to which the spin superstructure disap- 
pears somewhere between 35 mK and 1.5 K. Evidently, 
the transition from a rhomboid spin solid to the spin gas 
occurs at a temperature below 0.63 K, our lowest tem- 
perature for the plateau region. 

The solid line drawn in Fig. 2 through the 27.5 T data 
points below 5 K is the heat capacity of 1/8 mole of ideal 
two-dimensional Bose gas (3^ , with the particle mass be- 
ing the only fitting parameter. The fit yields a bosonic 
mass of (8.7±1.0) x 10~ 28 kg, which is 960 times the elec- 
tronic mass. It is 62% of the mass of the lowest triplet 



excitation, (1.4 ± 0.1) x 10~ 27 kg, obtained from a si- 
nusoidal fit of the dispersion measured by inelastic neu- 
tron scattering at zero field 33|. The similarity of the 
two masses suggests that the triplet excitations at lower 
fields are somewhat extended objects |27J similar to the 
ground-state triplets at the magnetization plateau. It 
is important to note that the thermal de Broglie wave- 
length of the spin gas at 0.63 K is as large as 1.7 times 
the average interparticle distance. Hence, the freezing of 
the spin gas into a rhomboid structure takes place in a 
highly degenerate regime. 

Surprisingly, the specific heat between 29.6 T and 33 T, 
in the field region above the magnetization plateau, is 
completely different from that in the region below the 
plateau and is identical with the 27.5 T data except for 
the low-temperature anomaly that signals phase tran- 
sition. In particular, it remains at R/8 in the mid- 
temperature range within our experimental accuracy, 
demonstrating that the same number of particles form a 
two-dimensional spin gas at these fields as at the plateau. 
It is natural to conclude that these particles at the two 
field regions are identical. 

On the other hand, the magnetization above the 
plateau increases continuously with increasing field and 
reaches a 52% higher value at 33 T 0, 0|. It is very 
unlikely that the extra magnetization is carried by the 
particles of the spin gas, since (S z ) will then have to be 
larger than unity and this would require mixing of an 
S z = +2 excited state. A presence of such a state should 
result in a low-temperature specific-heat hump whose po- 
sition changes with magnetic field, very much like the 
behavior shown in Fig. 1 for the field region below the 
magnetization plateau. We propose that the extra mag- 
netization is carried by additional S z ~ +1 particles that 
are thermodynamically inert at least up to the highest 
temperature of this experiment. 

Phenomenologically, the particle-number conservation 
that underlies the two-dimensional gas behavior can be 
explained by large energy required for excitations that do 
not conserve the particle number. We believe, however, 
that this remarkable phenomenon deserves a fundamen- 
tal explanation, which is presently lacking. 

Momoi and Totsuka [34| have predicted that two 
species of S z = +1 bosons compose the state above the 
1/3 magnetization plateau, beyond the field region ex- 
plored by experiments to date. They are those same par- 
ticles that form the plateau phase and additional par- 
ticles that condense into a magnetic superfluid, whose 
density increases with increasing magnetic field. A sim- 
ilar two-species state in the field region above the 1/8 
magnetization plateau can explain why the number of 
particles that make up the spin gas remains the same 
as in the plateau phase. However, the robustness of the 
particle number against thermal excitations still remains 
a mystery. 

The sharp transition marked by the low-temperature 
peak in the specific heat is probably due to freezing of the 
spin gas [3{| . The strong field dependence of the transi- 
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tion temperature indicates that either the particle mass 
or the strength of the repulsive interaction increases with 
the density of the thermodynamically invisible particles. 
Since freezing in two dimensions is generally a Kosterlitz- 
Thouless transition 36], which has no signature in spe- 
cific heat, the sharp peak suggests an importance of spin- 
phonon coupling or inter-layer spin coupling. 

In summary, we have uncovered a distinct spin be- 
havior in each of three high magnetic-field regions of 
SrCu 2 (B0 3 ) 2 . In the region below that of the 1/8 mag- 
netization plateau, where translational symmetry exists, 
the spin gap remains open, indicating hybridization of 
at least one of the two S z = +1 triplet states with the 
singlet ground state, which now acquires a magnetic mo- 
ment. We conclude that there is no phase transition be- 
tween this region and the low-field, non-magnetic region. 
The regions at and above the magnetization plateau con- 
stitute a second contiguous phase, where those triplets 
whose number equals 1/8 of the spin dimers form a 
two-dimensional gas of massive bosons above the field- 
dependent freezing temperature. In addition, the region 



above the plateau contains thermodynamically inert par- 
ticles that carry the magnetization in excess of 1 /8 of the 
saturation value. Further studies are called for to identify 
these carriers and to elucidate their role in the ordering. 
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